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Band-alignment induced current modulation in Bi2Se3 three-dimensional topological insulator slab has 
been investigated by quantum transport simulations for three different device designs, one for purely lateral 
transport and other two with vertical transport. Non-Equilibrium Green Function formalism has been 
deployed to understand the transport mechanism in band- alignment devices to appraise the possibility of a 
3D-TI based resonant device. A resonance condition is observed when the Dirac-points (bands) are aligned. 
This results in the maximum current at resonance for the design with only lateral transport. However, 
current ratio between resonant and non-resonant condition is found to be relatively small and strong 
temperature dependence is also noticed. The other two designs with vertical transport have degraded 
transfer characteristics, although from state-of-art literature they are expected to manifest nearly an ideal 
resonance peak. The physical insights for these observations have been posited along with the suggestions 
for attaining close to an ideal operation for the first design, which we also suggest for the pursuit in the future 
for spintronic oscillators and analog multipliers based on band- alignment induced resonance. 

Graphene-like gapless linear Dirac-bands, but in odd numbers, on the surface of three-dimensional (3D)^ 
topological insulators (TI), with bulk insulating phase, have excellent transport properties. Time-reversal 
symmetry protects these surface states against any non-magnetic perturbation^ to ensure high conduc- 
tivity by suppressing backscattering, which makes 3D-TI suitable for applications like interconnects^ and 
quantum computing^. Furthermore, spin-momentum locking^ for carrier transport in these helical surface states 
projects TI also as a promising candidate for spintronic applications^'^, which gives it an advantage over spin- 
degenerate dual-valley graphene Dirac-bands. Among the extant 3D-TI, BiaSea^ has the largest bulk bandgap of 
0.3 eV^ which provides large energy window for exploiting these exotic surface states in TI, and hence has been 
chosen as representative material in this work. Recently, there has been a surge in device proposals^°"^^ to regulate 
electron flow in Dirac-band materials via band-modulation by filtering out unwanted momentum states. 
Furthermore, by leveraging on the Dirac-bands of graphene, recently there has been a proposal of the symmetric 
Field Effect Transistor (SymFET)^^, a novel device with a resonant current peak for an appropriate balance of 
channel bias (Vds) with Fermi-level {E-f). This device based on vertical tunneling from n-doped layer to p-doped 
layer, separated by an insulator, is expected to exhibit a large current switching ratio which is symmetrical about 
resonance peak, at which the Dirac-points (bands) exactly align between two layers (see Fig. 4(b) of Ref. 13). The 
device operation has also been shown analytically^^ to be robust against temperature variation. Consequently, this 
device has been considered as a good candidate for designing high-speed analog multipliers and oscillators^l Such 
high-speed multipliers and oscillators in spintronic circuits would serve an important role in hybrid logic based 
systems^^ with 3D-TI serving a common platform for components. 

In this work, therefore, we comprehensively investigate the transport mechanism in band- alignment devices to 
appraise the possibility of a similar resonance state in 3D-TI which may serve as an analogue of multipliers and 
oscillators in the spintronic circuits. For this purpose three different device designs are evaluated, one purely for 
lateral transport (Design- A) and two for vertical transport (Design-B and Design-C), to operate in the helical 
surface bands. Design-A presents a simple understanding of band- alignment concept and explains the effect of 
mode- matching (electron momentum modes), gate-tunable density of states (DOS) and temperature on transfer 
characteristics of the device. Design-B emulates the symFET concept in 3D-TI to peruse its assumptions via 
quantum transport investigation. Although a similar resonant state is observed, its effectiveness is much weaker 
than originally expected for which reasons have been discussed. Design-C takes the ideal symFET operation to a 
more practical case on 3D-TI. However, there is a complete loss of resonance in Design-C. Therefore, best results 
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Figure 1 | (a) Device Structure (Design -A) for quantum transport modeling through Bi2Se3 slab with shorted top and bottom gates, (b) Potential 
distribution along the transport direction, and i^d are electrochemical potential at the source and drain end respectively. Gate voltage Vtg shifts the 
potential energy in the channel. A represents the energy difference between and the Dirac-point in the source contact and is used as a reference energy 
scale throughout this work, (c, d, e) Schematic of the band profile along the transport direction (source, channel and drain from left to right) for 
three different gate voltages, illustrating the mechanism of mode matching. At 0 K transport is within the energy window delineated with blues dot- 
dashed line. Blue arrows depict ballistic transport along the transport direction with conserved transverse (along y-axis) momentum. Green bands in (d) 
emphasize the band-aligned condition between source and channel region. 



are obtained for Design-A although the current ratio, between res- 
onance state (Ir) and non- resonance state at zero gate voltage (Iq), is 
still relatively small and the device transfer characteristic is sensitive 
to temperature. Accordingly, suggestions to further improve the ratio 
for Design-A are provided. 

The work is organized as follows. In methods section, the device 
Hamiltonian is introduced and the methodology for computing 
transport is explained. Results and Discussion section discusses the 
results of our simulations and unravels the effect of band-alignment 
operation on the carrier transport in 3D-TI devices. This is followed 
with the summary and conclusions of our study. The supplementary 
material provides the mathematical proofs of some important con- 
cepts discussed in this work for more rigorous readers. 

Methods 

For describing the device Hamiltonian, the parameters extracted from ab-initio cal- 
culations, are optimized for both in-layer and inter-layer coupling to fit the experi- 
mentally observed energy- dispersion for Bi2Se3 as elaborated in our previous work^^. 
Each infinite quintuple layer (QL ~ 0.943 nm)*^ (x-y plane; see Fig. 1(a)) of Bi2Se3 is 
described by basis set of up and down orbital of two equivalent Se atoms in k.p 
model by Hp as. 
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where v = 2.5 eV-A is the Fermi velocity, = 0.125 eV~^A~^ and m2 = 
— 0.04 eV~^A"^ are the orbital masses and the parameter d = —0.22 eV is introduced 
to generate a gap. The hopping between adjacent layers (z-axis) is described by tight- 
binding parameter = 0.35 eV in hopping matrix T as. 
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In the x-direction, kx wave-vector in Eq. (1) is discretized by substituting it with -i 
d/dx. Transverse direction (y-axis) is modeled in k-space (eigenmodes) assuming 
periodic condition (infinite width) because there is no potential variation along y- 
axis^^. Non-Equilibrium Green function (NEGF) algorithm^^ with Recursive- Green 
function^** is implemented to describe the ballistic transport through a defect-free 3D- 
TI device. Ideal contacts are modeled as semi-infinite channel region and the contact 
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self-energ)^ is computed from self- consistent solution of surface green function^^. The 
spatial distribution of inhomogeneous potential in the all three devices, because of 
concurrent application of both gate and channel bias and the potential drop at contact 
leads, has also been modeled within NEGF formalism. The detailed potential distri- 
bution for Design-A and Design-B-C has been elaborated in results and discussion 
section for lateral and vertical transport, respectively. Next, the gates are phenom- 
enologically modeled with a translation factor of one to regulate the potential in the 
channel. We note that self- consistent solution of channel potential with a Poisson 
solver would disturb the perfect band-aligned resonant condition. Since, this would 
both degrade the results and eclipse the transport concepts which we in particular 
want to stress upon, the Poisson equation has not been exercised in this work. 
However, if the device performance is to be examined for practical implementation, 
the NEGF transport must be solved self- consistently with Poisson to correct the 
charge induced potential variation in the system. 

Results and Discussion 

Transport behaviour in three different band- alignment device 
designs is expounded in this section. For Design-A, the transport is 
only in lateral direction which makes it easy to understand the band- 
alignment concept. It is furthermore supported with derivations in 
supplementary section SI and S2 to provide clear insights into the 
mechanism. Next we delve into more complex S-shape devices for 
vertical transport to be in-line with a similar symFET^^'^^ proposal on 
graphene. The transport mechanism and limitations have also been 
discussed for each design. 

Band- Alignment for Lateral Transport. Figure 1(a) illustrates the 
device structure for Design-A with 13 QL thickness (z-axis) and 
30 nm channel length (x-axis) along with semi-infmite contact on 
either sides of the channel. The sample thickness has been chosen to 
be sufficiently large so that inter- surface coupling effects are 
negligible^ ^ because the surface wavefunction decays to nearly zero 
beyond 10 QL^^ Top and bottom gates are shorted to have uniform 
influence on all layers of 3D-TI slab. The potential profile along the 
transport direction is shown in Fig. 1(b) to illustrate the equal drop of 
bias voltage (Vds) at the source and drain terminals and a flat 
potential in the channel, which shifts the energy-bands in the 
channel and drain region by -Vds/2 and -Vds respectively. An 
important conceptual point to note here is that the energy bands 
are strictly defined only in the equilibrium condition i.e. in the 
absence of channel bias. However, in non-equilbrium, under- 
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Figure 2 | Design-A. (a) Operation for 13 QL thick 3D-TI for three different Fermi-Levels in the surface bands (Vds = 2 A) at 0 K. Observe the 
resonant peak at — A. Ir/Iq ratio is defined between Vtg = 0 (lo) and Vxg = — A (Ir) and stated over each trendline with respective color, 
(b-d) Transmission spectrum (in per [im) for Vtg equal to 0 (b), —A (c) and — 2A (d) at 0 K for A = 0.04 eV as a representative case from (a). 



standing in terms of the shift of energy-bands along the transport 
direction to follow quasi Fermi-level distribution ia a good appro- 
ximation and provides simple insights into transport behaviour (see 
M. Lundstrom, Fundamentals of Carrier Transport: Cambridge 
University Press, 2009). Next, also note that because of selenide 
vacancies^^ a typical sample of Bi2Se3 TI is n- doped and a small 
bias applied to electrically characterize the material for ballistic 
transport should approximately follow potential distribution of 
Fig. 1(b). Application of gate voltage over such a system would 
result in a commonly observed ambipolar V-shaped transfer charac- 
teristics^^ which are just a signature of Dirac-bands. However, despite 
of similar device geometry the band- alignment device operates 
differently as shown in this work. 

To assay the transfer characteristics of band- alignment device, a 
channel bias of V^s — 2 A (see Fig. 1) has been applied across the 
device, where A represents the energy difference between [is ^nd the 
Dirac-point in source contact and is used as a reference energy scale 
throughout this work. Therefore, in equilibrium condition A equals 
Fermi-level (£y), with zero-energy level calibrated at the Dirac-point. 
Energy-bands in the channel are shifted by applying the gate voltage 
(Vtg)> which configures it from a purely n- doped (Vtg ^ 0 volts) to 
purely p-doped (Vtg = — 2A volts) state (see Fig. l(c-e)). The res- 
onant condition in this device architecture is achieved at Vtg = ~ A, 
when the bands in the source region are exactly aligned with the 
bands in the channel so that wavevector modes along both x-axis 
(kx) and y-axis (k^) are conserved resulting in perfect transmission 
across the source- channel interface. Also, observe that there is a 
mode-mismatch across the channel-drain interface, and degrades 
as the Vtg is swept from zero to —2 A. This effect manifests clearly 
in the transport results which we discuss next. 

Figure 2(a) illustrates the transfer characteristics of Design-A. It 
shows a clear evidence of resonant current peaks at Vtg = ~ A (see 
Fig. 1(d)), for different Fermi-levels (Ef). For operation in the linear 
Dirac-bands, the current increases with the increasing magnitude of 
Ef because the density of states (DOS) scales linearly with £y ^ for 
Dirac bands, which subsequently increases the current (see Eq. 
(Sl.l)). Furthermore, Ir/Iq ratio which is defined between Vtg — 0 
and Vtg = ~A also increases with the current and becomes more 
clear at higher Ef. However, note that since resonance condition 
requires Vtg = ~ A and Vds — 2 A, further increase in £y will drive 
transport energy range into the bulk bands, while for spintronic 



applications we need to operate in the surface-bands to leverage 
the spin-momentum locking for electrons. This also shows that a 
new 3D-TI material in future with a larger bulk bandgap may 
improve the ratio. The transport characteristics are further illustrated 
via transmission spectrums in Fig. 2(b-d). Specifically notice the 
band- alignment with nearly symmetric operation about the Dirac- 
point for Vtg = — A in Fig. 2(c). Next, since the transverse modes (ky 
modes) for carrier propagation in the device, which is uniform along 
the y-axis, should always be conserved^^, the transit from the larger to 
smaller modes is prohibited and they are filtered out (see Fig. 2 of Ref 
26). Therefore, light bluish-yellow region in Fig. 2(b) and 2(d) show 
significant mode filtering in valence and conduction band region 
respectively for Vtg equal to 0 and —2 A. Mode-filtering also clearly 
manifests for valence band region in Fig. 2(c) and 2(d) because of the 
mode-mismatch between channel and drain (see Fig. l(c-e)). 
Although, the mode-filtering may seem unlikely for resonance con- 
dition in Fig. 2(c), note that unlike graphene, the conduction and 
valence bands of Bi2Se3 3D-TI are not symmetric^'^^ Valence band 
has smaller Fermi-velocity than conduction band. Therefore, there is 
mode-filtering in low-energy range even at resonance because elec- 
trons attempt to transit from large ky modes of valence band in the 
channel to conduction band in the drain (see Fig. 1(d)) which does 
not support those modes. There are two additional observations 
from the transfer characteristics in Fig. 2(a). Firstly, the current at 
two extremes of Vtg (0 and —2 A) is not equal (higher at Vtg = 0)- 
Secondly, there is a local minimum on either sides of resonance. The 
first happens because of mode-mismatch (see supplementary SI) 
between bands in the channel and drain and therefore is strongest 
for Vtg — —2 A. The effective DOS in the transport energy window 
between |is and jio (see Fig. 1) is responsible for the second obser- 
vation, which is mathematically proven in supplementary S2. As the 
Vtg is swept from — A to 0 (or —2 A), the potential in the channel 
starts to move the higher DOS out of energy window on one end of 
the energy-grid and replaces it with lower DOS at another end of the 
grid, which should have resulted in continuous decrease in current. 
However, the effect of DOS competes with the mode matching 
between the source and channel and results in actually higher current 
at Vtg = ~A. The inflection points mark the transition of dom- 
inance of one effect over another. 

The effect of temperature on operation of Design-A is next 
appraised in Fig. 3. The temperature results in the spread of 
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Figure 3 | Design-A. (a) Current vs gate voltage at 0 K (blue circles), 100 K (red square), 200 K (magenta triangle) and 300 K (black stars) for A = 
0.04 eV. Ir/Iq ratio and temperature are stated with respective colors for each trendline. Observe that signature of mode matching is severely weakened 
due to spread in Fermi- distribution, (b, c) Transmission spectrum (in per |am) for two gate voltages at 300 K. Compare it with Fig. 2 to observe the 
extinction of signature of mode match/mismatch in quantum transport. Dark red region is bulk conduction band. 



Fermi- distribution and therefore results in the inclusion of wider 
energy range and ky modes for the transport. This results in more 
number of energy and transverse modes which suffer mode-mis- 
match either at the source-channel or at the channel-drain interface, 
mainly the later one, which results in near complete loss of resonance 
at 300 K as illustrated in Fig. 3(a) by plotting the transfer character- 
istics and further confirmed via transmission spectrum in Fig. 3(b) 
and 3(c) respectively for Vtg equal to zero and — A. In the spectrums, 
observe that the two zero transmission points (blue dots) at zero kya 
correspond to the Dirac-points for bands in the source and drain 
terminal. The resonance condition for Vtg — ~^ can also be 
observed by comparing the width of bluish-green region between 
the two Dirac-points. Furthermore, the bulk conduction band (in 
red) is more clearly visible than the valence bulk bands (light green- 
yellow lines) because the Fermi-level is in the conduction band and 
since the contribution of energy points is weighted by the Fermi- 
distribution, the transmission is stronger in the bulk conduction 
band which are nearer to source electro -chemical potential than in 
the valence bands. Therefore, precisely the averaging of results over 
more mode-filtered wavevectors wanes the resonance at higher tem- 
peratures. Since the Ir/Iq ratio falls short of values that can be com- 
fortably measured in experiments, for instance around ten^^, it is 
important to enhance the ratio. To improve the Ir/Iq ratio a mode 
selective switch should be implemented in the channel or at the 
source-channel interface to completely filter out the non-normal 
modes (ky 0). This would in general force the device to operate 
at only one energy point between |is and i^d at which ky would be 



Design B and Design C 



conserved, but at resonance condition electrons would be able to 
conduct at all the energy points because the band-alignment would 
enable the ky conservation at all of the electron energies. This would 
result in large Ir/Iq ratio as expected for Graphene symFET in Ref 13. 
To achieve the suppression of non-normal modes, tilted gates or an 
insulating physical barrier like a lens can be implemented in the 
channeP^. Heterostructures^^ is another possibility which can be 
investigated in the future work. 

Band- Alignment for Vertical Transport. This section investigates 
two S-shape devices whose architecture is illustrated in Fig. 4(a). The 
semi-infinite source contact is coupled to 13 bottom QLs (bottom 
slab), while semi-infinite drain to top 13 QLs (top slab) with 4 QLs in 
the middle (middle slab). Large thickness of bottom and top slab 
ensure zero value of surface wavefunction in the middle region. 
Therefore, for operation solely in the surface bands, electrons must 
tunnel through the bulk bandgap to transit from source to drain. 
Next, for Design-B, the top and bottom gate voltages uniformly 
modulate the doping electronically in respective slabs. For the 
channel bias V^sj all layers of bottom (top) slab are hold at |is (iid) 
electro- chemical potential of the source (drain), while a linear 
gradient (vertically) is considered in the middle. Note that the 
physical isolation of layers with contacts on either ends should 
hold the potential closer to respective contact electrochemical 
potentials^°'^\ in contrast to Design-A in which both sides of the 
channel were equally coupled to contacts across all layers and 
hence had equal drop at both ends (see Fig. 1(b)). Nonetheless, the 
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Figure 4 | (a) Device Structure for vertical transport in 3D-TI. Source is coupled to bottom few layers (bottom slab) and drain to top few layers (top slab) . 
As a representative case, bottom and top slabs are 13 QL each, whereas middle region is 4 QL. Electrons injected from the source have to tunnel through 
the middle slab for operation in the bulk band gap. For Design-B, the top and bottom gate voltages uniformly modulate the doping electronically 
in the respective slabs. For Design-C, a more practical scenario for 3D-TI is presented with uniformly n-doped sample, (b) Potential distribution in x-z 
plane for band-aligned condition. Blue arrow depicts vertical transport from Dirac-bands of bottom slab to that of top slab, tunneling through the middle 
region. For channel bias, all layers of bottom (top) slab are hold at Hs (|^d) electro-chemical potential of source (drain), while a linear gradient is 
considered in the middle. Gradient for gate potential in Design-B is assumed to be only in the middle region, whereas for Design-C it is from bottom to top 
layer of the device (across all three regions). 
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Figure 5 | Design-B. (a) Current vs top gate voltage at 0 K for A = 0.02 eV (red squares) and A = 0.04 eV (blue circles). Observe the resonant 
condition at Vjc — —2 A. (b) Illustration of current flow in the device. Observe that the current mainly flows laterally (along x-axis) on top and bottom 
surface but is forced to transit vertically (z-axis) at the device boundary, (c-h) Transmission spectrum (in per [im) for Vjc equal to 0 (a), — A (b), — 1.5A, 
(c) -2A (d), -3A (e) and -4A (f) at 0 K for A = 0.04 eV as a representative case from (a). 



precise drop at contacts should be computed by solving poisson 
equation self- consistently which has been explicitly not done as 
explained in methods section. Gradient for gate potential in 
Design-B is assumed to be only in the middle region which 
ensures degeneracy for both the Dirac-bands in source and drain 
region each. The potential-profile for Design-B warrants that there 
is absolutely no mode mismatch between bottom (top) layer surface 
band of source (drain) and bottom (top) slab of the channel. The 
bottom gate is grounded and top gate voltage is applied to achieve 
band- alignment condition of resonance which is now expected at 
Vds — 2 A and Vtg = —2A because potential in top and bottom 
slab is assumed to be on par with electrochemical potential of the 
contacts, instead of fifty-percent drop as considered for Design-A. 
Therefore, for operation of Design-B, doping in the top slab, 
including drain contact to eliminate mode-mismatch at the top 
slab-drain interface, is regulated by an electronic equivalent of Vtg 
value (see Fig. 4(b)). This profile is necessary to emulate the scenario 
similar to graphene symFET^° for perfect resonant transmission 
from the bottom Dirac-band to the top Dirac-band. For Design-C, 
a more practical scenario for 3D-TI is presented with a uniformly n- 
doped (de-facto sample) material because experimentally, at least at 
present to the best of our knowledge, it should be difficult to precisely 
control the doping to specific layers in 3D-TI. The gate potential 
gradient is taken from bottom to top layer of device (across all 
three regions) (see Fig. 4(b)) to follow an expected setup. Such a 
profile, however introduces mode-mismatch at the top slab-drain 
interface, as well as, at the top layer of the source contact and the 
bottom slab. The purpose of having two device designs is to first via 
Design-B examine the symFET operation comprehensively through 
quantum transport, to go beyond the analytical evaluation in Refs. 
13, 20, and then via Design-C assess its feasibility in 3D-TI because 
the relaxed and more realistic gate potential gradient in an n-doped 
sample may deteriorate band-alignment operation. 

Figure 5(a) illustrates the transfer characteristics of Design-B. It 
shows a clear evidence of resonant current peaks at Vtg — —2 A, for 
two different A values. Since, larger A also enforces higher Fermi- 
level for the band-alignment devices discussed in this work, a larger 
current is observed for A = 0.04 eV (red-squares) than for A = 
0.02 eV (blue circles) because higher Ef also results in transport 
through energy range with larger DOS (see Eq. (15) of Ref 24 or 
Eq. (S2.1)). However, interestingly Fig. 5(a) also illustrates very 
strong asymmetry about resonant condition, with larger current 



for more electronegative gate potentials in strong contrast to 
Design-A which shows lower current for similar condition. To ana- 
lyze this anomaly. Fig. 5(b) illustrates the current spectrum in the 
device. Its ordinate axis represents the QLs of the device (z-axis) 
while abscissa denotes lateral direction (x-axis). The reddish-yellow 
color shows the region where current is mainly concentrated in the 
device while arrows mark the region where most current transits 
vertically from bottom to top slab. Observe that current flows ver- 
tically mainly at the device boundary, otherwise it flows laterally (in- 
plane) along X-axis. This is in sharp contrast to graphene symFET in 
which the transport is assumed to be only along the z-axis^^ In 
symFET, if the current flow is only along the z-axis, then both kx 
and ky wavevectors together constitute the transverse modes (net 
magnitude equals root of sum of square of both wavevectors), and 
their conservations acquiesce just one energy point between |is and 
jiD for conduction in the non- resonant state and through the entire 
energy grid for resonant condition resulting in a sharp resonance 
peak. Since, NEGF quantum transport simulation in Fig. 5(b) shows 
current flow along x-axis as well, the result implies that only ky 
wavevector must be conserved, analogous to Design-A. A strong 
resonance peak therefore may not be expected for such a system. 
Other implications from current distribution in Fig. 5(b) are that 
firstly the transfer characteristics should strongly be regulated by 
lateral transport and secondly the vertical tunneling based band- 
alignment effect may be weaker than initially anticipated. Never- 
theless, note that the perfect mode-matching atleast at the contact- 
channel interface for both the top and bottom slabs simplifies the 
transport mechanism with respect to Design-A. Furthermore, to 
advance the understanding of the asymmetry in transfer character- 
istics. Transmission spectrums for a few Vtg values are illustrated 
from Fig. 5(c-h) for increasingly negative values. The dark blue spot 
at 0 eV energy in all spectrums corresponds to the Dirac band for the 
source contact and the bottom slab, whereas the other dark blue dot 
for ky = 0 in Fig. 5(d-h) is consistent with movement of the Dirac- 
band in the drain contact and the top slab (see Fig. 4(b)). For analyz- 
ing these spectrums, recall that valence bands in 3D-TI have lower 
Fermi-velocity than in conduction bands. Therefore, DOS is higher 
in valence band region (see Eq. (S2.1)), which results in wider energy- 
dispersion. For Fig. 5(c), the transport in low energy range is from 
wider valence band region to more steep conduction band region, 
which results in mode-filtering (see clipped light blue region in low 
energy range). All transverse modes are perfectly filtered out at the 
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Figure 6 | Design-C. (a) Current vs top gate voltage at 0 K for A = 0.02 eV (red squares) and A = 0.04 eV (blue circles). Observe the absence of 
resonance condition at Vjc — —2 A. (b-f) Transmission spectrum (in per [im) for Vjq equal to 0 (b), — A (c), — 2A (d), — 3A (e) and — 4A (f) at 0 K for 
A = 0.04 eV as a representative case from (a). 



energy value corresponding to Dirac-point in the top-slab. Red 
region in spectrum for Vtg = — 2A in Fig. 5(f) shows the band- 
alignment condition which clearly evinces the largest transmission 
in comparison to Fig. 5(d), 5(e) and 5(g). However, also recall from 
the discussion for Design- A and the proof in supplementary SI that a 
band-aligned condition results in the lowest DOS configuration in 
transport energy window. The lateral transport in top layer in high 
DOS region of conduction band for Vtg ^ 0 and even higher DOS of 
valence band for Vtg — —4 A, results in even higher transmission 
(see Fig. 5(c) and 5(h)) and current than the resonant condition. 
Larger DOS in valence bands is also a precursor for the asymmetry 
observed in Fig. 5(a) (see Eq. (Sl.l)). Therefore, although we observe 
resonant condition in Design-B, the multi-mode conduction and 
dominance of DOS effect due to lateral transport does not position 
Design-B as a switch. 

Figure 6(a) illustrates the transfer characteristics of Design-C for 
two different A values and does not show any resonance behavior. 
The current is higher for larger A for the same reason as Design-B 
above. However, the lowest current point is not at Vtg ^ — 2A as 
would have been expected for operation in low DOS configuration 
for the band-aligned condition. The current distribution in Design-C 
is qualitatively same as illustrated for Design-B in Fig. 5(b); vertical 
transport is chiefly at the device boundary. Therefore we follow a 
similar approach to understand these observations, by investigating 
DOS and mode-matching via Transmission spectrum in Fig. 6 (b-f) 
for increasingly electronegative Vtg- For Fig- 6(b) at Vtg ^ 0» there is 
absolutely no difference between Design-B and Design-C, however 
subsequent application of Vtg in Design-C does not affect the dop- 
ing in drain contact and it retains same doping as the source terminal. 
Therefore, Vtg now introduces strong mode-mismatch at the drain- 
channel interface which degrades further for the larger magnitudes of 
Vtg» a problem similar to Design-A. Besides, the DOS in the trans- 
port energy window and band- alignment for Dirac-bands of the top 
and bottom slab vie for the effective results. For Vtg — ~ A, bands in 
the bottom slab, top slab and drain are positioned similar to Fig. 1(c) 
for Design-A and therefore Fig. 6(c) shows mode-filtering through 
yellowish arc in the valence band similar to Fig. 2(b). The surround- 
ing light blue region results from near-perfect lateral transport from 
the source contact to the bottom slab because there is no mode- 
mismatch at bottom most layer of the device, although there is a 



minor mismatch at top of the bottom slab due to the gate-potential 
gradient. Next, Fig. 6(d) illustrates the band- alignment condition in 
Design-C (see Fig. 1(d)). It is this band- alignment operation between 
the bottom and top slab which results in slightly higher current than 
Vtg — —1.5 A in Fig. 6(a), for an otherwise operation in the lowest 
DOS configuration (see note for Design-C in supplementary S2). 
Fig. 6(e) illustrates mode-filtering scenario in Design-C, a case sim- 
ilar to Fig. 1(e) and its corresponding spectrum in Fig. 2(d) of Design- 
A (compare yellow arc in the conduction band). The light blue region 
is again due to near-perfect lateral transport in the bottom slab. As 
Vtg is further increased, the red-region in Fig. 6(f) reveals the high 
DOS states of valence band and predicts an increase in the current 
through the device which is observed in Fig. 6(a). Therefore, from 
transmission spectrum we understand that Design-C inherits few 
characteristics of both Design-A and Design-B, in which vertical 
and lateral transport, band-alignment between the top and bottom 
slab and mode-mismatch at contact- channel interface have to be 
analyzed simultaneously to understand its transfer characteristics. 

Finally, we would like to note that our model indeed shows side 
surfaces as illustrated in Fig. 1 of our previous work^^. However, since 
the thickness of each slab, we consider, is only 13 QL (and additional 
intermediate slab of only 4 QL for Design-B and C), the strong 
confinement effect along the z-axis destroys topological side surface 
states. The previous studies on side-surface for energy- dispersion^^ 
and transport^^ had considered this axis to be infinite or semi-infi- 
nite. Similarly, experimental works to claim evidence of quantum 
transport through side- surfaces were performed on extremely thick 
Bi2Se3 nanowire (140 nm in Ref 34 and 200 nm in Ref 35). In fact, if 
a similar strong confinement were applied along the y-axis, the topo- 
logical states on the top and bottom surfaces should also be 
destroyed. 

Conclusion 

In summary, this work expatiates on the concept of band- alignment 
or mode-matching induced resonance to examine the feasibility of a 
resonant device for Bi2Se3 3D-TI, a spintronic material. The three 
different device-designs are evaluated for this purpose. Design-A, 
which is based only on lateral transport, shows that operation in 
low DOS energy-range at resonance competes with mode-matching 
and mode-filtering to manifest the observed transfer characteristics. 
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It is found that despite of a resonance peak, the Ir/Iq ratio is not very 
exciting. The degradation of the ratio is chiefly because the transport 
in non-resonant condition is not exactly Umited to one energy point 
between contact electro -chemical potentials (jis and iid) as may be 
expected for Graphene symFET. Therefore, the solutions to limit the 
transverse modes for improving the device performance have been 
suggested for future improvements. It is also observed that the rise in 
temperature degrades the resonance behaviour because of averaging 
over wider energy and ky modes with mode-mismatch. In Design-B 
and Design-C the vertical and lateral transport, band- alignment 
between the top and bottom slab and mode-mismatch at contact- 
channel interface have to be analyzed simultaneously to understand 
its transfer characteristics. Design-B although exhibits a resonance 
peak, is found to be less effective than Design-A. Design-C which 
could be a more realistic implementation on 3D-TI; there is a com- 
plete loss of resonance peak. Therefore, as a future direction we 
suggest to further improve the Design-A design to enhance the ratio 
for implementing spintronic circuits using 3D-TI. 
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